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摘要 :基于 2000 一 2014 年 的 MODIS-NDVI 与 MODIS-LST 数据 ,利用 温度 植被 干旱 指数 对 黄土 高 原 土壤 湿度 进行 了 反 演 ,采用 时 
滞 互 相关 法 分 析 了 土壤 湿度 与 植被 覆盖 的 年 内 变化 特征 及 其 时 空 相互 关系 。 结 果 如 下 :(1) 黄土 高 原 植被 生长 对 土壤 湿度 变 
化 存在 明显 的 时 请 效应 ,植被 生长 以 及 植被 物候 与 土壤 湿度 变化 密切 相关 。(2) 黄土 高 原 土壤 湿度 对 其 植被 覆盖 变化 影响 强 
烈 区 ,主要 分 布 在 农业 与 草原 生态 区 北部 的 丘陵 沟 密 区 ,植被 生长 对 土壤 湿度 响应 迟缓 ;而 植被 覆盖 对 土壤 湿度 变化 影响 强烈 
区 ,分 布 于 西部 的 高 寒 地 带 ,响应 时 间 相 对 最 短 ,其 次 为 农业 与 草原 生态 区 内 河流 中 上 游 流 域 ,上 且 土 壤 湿度 对 植被 覆盖 的 啊 应 较 
快 , 自 东 向 西 . 自 北向 南 响应 逐渐 加 快 。(3) 从 植被 类 型 来 看 , I 类 植被 的 植被 覆盖 与 绰 坏 湿度 负 相 关 性 较 强 ,植被 生长 对 土壤 
湿度 变化 的 响应 程度 依次 减弱 ,响应 速度 较 慢 且 依 次 加 快 。 工 类 植被 的 植被 覆盖 与 土壤 湿度 的 相关 性 很 弱 ,滞后 时 间接 近 于 0 
天 ,反映 出 植被 覆盖 与 土壤 湿度 同步 变化 。 土 壤 湿度 变化 滞后 于 亚 类 植被 ( 旱 生 / 沙 生 植被 ) 的 植被 覆盖 变化 ,植被 覆盖 正 向 影 
响 土壤 湿度 变化 的 强度 依次 增加 ,但 时 间 效 应 依次 延长 。V 类 植被 (高 出 /高 寒 植 被 ) 生 长 对 土壤 湿度 变化 的 正 向 影响 程度 最 
强 , 而 且 响 应 速度 较 快 。(4) 黄土 高 原 植被 生长 与 土壤 水 分 的 年 内 变化 具有 密切 的 关系 ,土壤 水 分 年 内 变化 的 峰值 和 谷 值 与 植 
被 物候 期 非常 吻合 ,因此 土壤 水 分 的 年 内 变化 可 作为 植被 物候 特征 提取 的 一 种 重要 依据 。 
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Abstract: Based on the MODIS-NDVI and MODIS-LST data from 2001 to 2014, we have estimated the soil moisture using 
temperature vegetation. dryness index, analyzed the spatiotemporal seasonal variation of the vegetation cover and soil 
moisture and their\interrelationship in the Loess Plateau. The findings suggest the following: (1) There are obvious time-lag 
effect for vegetation growth relative to soil moisture in the Loess Plateau. The vegetation cover and vegetation phenology is 
closely-related to soil moisture change. (2) The intense affected area of soil moisture to vegetation cover located mainly in 
the mountains and hills of the northern agriculture and grassland biome, which the vegetation cover has slow response. The 
intense affected area of vegetation cover to soil moisture spread over the western alpine region, which the response time of 
soll nioisture ls relatively shortest. The second one distributes across the upper and middle reaches of the rivers in the 
northern agriculture and grassland biome, and the soil moisture respond slowly. The response speed gradually increase from 


east to west and from north to south. (3) From vegetation types, there are relatively strong negative correlations between 
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vegetation cover and soil moisture of vegetation | that the response degree of vegetation growth to soil moisture change 
weakens successively, and their response rate ls relatively slow and speed up. A weak correlations were found between 
vegetation cover and soil moisture of vegetation lH , and their lag time are close to 0 day, which reflects vegetation cover and 
soil moisture change simultaneously. The soil moisture changes lag behind vegetation cover change of vegetation Ll 
( xerophytic/psammophytic vegetation ) , and the degree that vegetation cover positive impacts on soil moisture increase 
consequently, but time effects extend. The positive impact of alpine vegetation ( vegetation lV ) to soil moisture is strong 
which has better response speed. (4) There is close relationship between the annual variance of vegetation growth and soil 
moisture in the Loess Plateau. The vegetation phenology ls a good fit to the peak and valley of soil moisture changes over the 


years which can be used as an important evidence for the extraction of vegetation phenological information . 


Key Words: soil moisture; vegetation cover; Whittaker smooth; time lag correlation; the Loess Plateau 


植被 覆盖 变化 是 区 域 生态 系统 环境 变化 的 重要 指示 ,是 衡量 地 表 植 被 状况 的 重要 指标 ;对 水 文 . 生 态 全 
球 变 化 千 都 具有 重要 意义 ;同时 ,植被 履 盖 又 是 有 影响 土壤 侵蚀 与 水 土 流失 的 主要 因子 汪 随 着 气温 的 不 断 升 高 ， 
地 表 蒸 散发 加 大 、 河 湖 萎缩 与 干 润 、 雪 线 抬升 ,气候 和 环境 的 干旱 化 成 为 严峻 的 生存 环境 问题 ,进而 使 得 土壤 
湿度 继续 降低 ,土壤 旱 化 程度 加 重 , 不 仅 反馈 给 植被 发 生 水 分 胁迫 ,致使 植被 正常 生理 代谢 失调 ,从 而 影响 地 
表 植 被 覆盖 的 变化 ,而 且 当 土壤 形成 不 可 逆转 的 干 层 时 , 叉 会 影响 到 降水 或 地 表 水 入 渗 补 给 地 下 水 ,以 及 陆地 
水 循环 ”。 由 此 可 见 ,土壤 水 分 是 气候 .土壤 和 植被 影响 水 分 循环 和 水 分 循环 影响 植被 生长 综合 起 来 的 关 
键 因 了 于 。 

土壤 水 分 是 植物 蒸腾 作用 和 光合 作用 的 限制 性 因子 ,土壤 水 分 决定 了 生态 建设 中 植被 类 型 和 结构 ,是 植 
物 生长 .植被 恢复 .土壤 侵蚀 的 主要 影响 因素 “而 黄土 高 原 地 处 干旱 一 半 干 旱 气 候 区 ,雨量 少 ,水 资源 缺 
乏 ,水 土 流失 和 殉 江 化 严重 ,生态 环境 十 分 脆弱 , 吻 破 坏 难以 恢复 ,加 之 土 层 厚 ,地 下 水 埋藏 深 ,土壤 水 分 是 植 
物 生长 所 需 水 分 的 主要 来 源 , 降 水 地表 水 及 地 下 水 转换 为 土壤 水 分 才 可 被 植被 吸收 ,因而 土壤 含水 量 的 多 窒 
与 植被 的 生长 密切 相关 ,进而 影响 植物 的 多 样 性 ; 义 通 过 影响 其 他 生境 因子 间接 地 影响 着 多 样 性 的 维持 ,而 生 
态 系 统 的 主要 “构件 ”物种 多 样 性 对 水 分 变 北 的 响应 直接 影响 着 生态 系统 的 功能 和 稳定 性 “" 。 因 此 ,土壤 水 
分 是 影响 黄土 高 原 环境 演变 的 重要 因子 ,作为 水 量 平衡 水循环、 生态 系统 演变 以 及 土壤 资源 的 重要 因素 ,已 
经 得 到 了 广泛 的 关注 “ 。 同 时 ,植被 覆盖 大 小 也 会 影响 土壤 湿度 发 生变 化 ,如 干旱 与 半 干 旱地 区 植被 覆盖 越 
大 根系 越 深 ,土壤 干燥 深度 与 程度 越 大 “” ,植被 覆盖 度 对 土壤 水 分 含量 影响 较 大 ,在 保持 其 原 有 的 丛生 密 根 
性 植物 建 群 和 较 高 覆盖 度 时 , 随 盖 度 增加 ,土壤 水 分 含量 呈 显 著 抛 物 线 型 趋势 增加 “” 。 刘 刚 等 研究 表明 除 农 
地 外 ,其 他 植被 类 型 下 都 不 同 程度 的 存在 着 土壤 干 层 , 干 层 严重 程度 ,林地 大 于 灌木 林地 ,灌木 林地 大 于 和 草地; 
草地 存在 临时 十 层 和 永久 干 屋 ”。 

土壤 湿度 的 传统 监测 方法 包括 中 子 仪 法 .快速 烘 干 法 .重量 法 .y 射线 衰减 法 .电阻 法 和 张力 计 法 等 ,这 些 
方法 能 有 效 监 测 圭 壤 湿 度 ,是 土壤 湿度 监测 研究 .应 用 的 主要 手段 ,但 是 该 方法 采样 速度 慢 .耗费 大 量 人 力 和 
物力 .使 用 范围 有 限 ,加 之 测 站 分 布 不 均 ,难以 实现 大 范围 实时、 精确 的 监测 。 随 着 遥感 技术 的 迅速 发 展 , 遥 
感 数据 的 多 时 相 、 多 光谱 ,多 分 辨 率 、 高 光谱 等 特性 ,使 得 土壤 水 分 快速 实时 动态 监测 成 为 可 能 “” 。 在 土 
壤 水 分 遥感 反 演 中 ,通过 地 表 温 度 和 植被 覆盖 特征 空间 反映 土壤 水 分 的 遥感 监测 得 到 了 广泛 应 用 ,温度 
植被 干旱 指数 (Temperature Vegetation Dryness Index, TVDI) 由 于 其 计算 过 程 输 入 数据 少 、 人 简单 易 行 .生物 物理 
学 含义 明确 等 特点 ,因此 本 文选 用 其 来 表征 土壤 水 分 的 高 低 。 本 研究 以 广义 的 黄土 高 原 范 围 作为 研究 对 象 ， 
基于 Matlab 软件 图 形 用 户 界面 (Graphical User Interface ,GUI) 功能 ,分 别 创建 了 TVDI 批量 生成 GUI、Whittaker 
平滑 GUI 和 时 沛 互相 关 GUI ,利用 植被 覆盖 与 地 表 温 度数 据 构 建 了 TVDI, 用 以 反映 黄土 高 原 土壤 湿度 ,并 分 
析 了 植被 黎 盖 与 土壤 湿度 的 年 内 变化 及 其 相互 关系 时 空 变 化 特征 ,进而 统计 不 同 植被 类 型 植被 履 盖 与 土壤 湿 
度 的 年 内 相互 关系 。 
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1 研究 区 概况 


对 于 黄土 高 原 的 范围 界定 问题 ,一 般 具 有 狭义 .中 
义 和 广 义 之 分 。 黄 土 高 原 的 狭义 范围 指 的 是 长 城 沿线 
以 南 ,秦岭 以 北 , 吕 梁山 以 西 ,青藏 高 原 东 缘 ,测算 大 致 
30 万 km ,强调 完整 的 自然 单元 ;中 义 范 围 指 长 城 沿 线 
以 南 ,秦岭 以 北 ,太行 出 以 西 , 日 月 山 以 东 ( 包 括 日 月 山 
和 青藏 高 原 东 缘 之 间 的 区 域 ) ,测算 大 致 43 万 km , 强 
调 的 是 完整 的 地 貌 单 元 ;黄土 高 原 的 广义 范围 为 阴山 以 
南 西北 以 贺兰山 为 界线 ,秦岭 以 北 , 以 东 包括 太行 山东 
态 , 日 月 山 以 东 , 强 调整 个 黄土 高 原 与 黄河 中 游 整 个 水 
土 流失 的 联系 。 


ee 
1000 2000 km 
| | 


本 文采 用 的 是 广义 的 黄土 高 原 范围 ,广义 黄土 高 原 | 
西 起 祁连山 支脉 马 峭 岭 东 迄 太行 山 , 北 到 长 城 y 南 抵 秦 Fig.1 The survey of research area 


岭 。 横 路 青 、 甘 、 宁 、 内 蒙 主 、 陕 . 晋 . 驳 7 省 区 大 部 或 一 

部 ,面积 约 63 万 km , 约 占 全 国 陆地 总 面积 的 8%。 是 东南 湿润 季风 气候 向 西北 内 陆 干 旱 气 修 过 渡 暧 温 沉 落 
阔 叶 林 向 典型 草原 和 荒漠 草原 过 渡 的 过 渡 地 带 。 受 温带 大 陆 性 气候 的 控制 ,气温 ,降水 量 季 相 分 明 ; 从 东南 向 
西北 ,气候 依次 为 半 湿 润 气 候 、 半 干旱 气候 和 干旱 气 修 ;土壤 依次 为 宰 士 \ 坊 二 \ 黄 绵 土 和 灰 钙 土 ;植被 包含 11 
个 植被 型 组 ,23 个 植被 型 ,171 个 群 系 ,11 个 植被 型 组 分 别 为 针 吐 林 < 沙 叶 林 、 灌 从 、 草 从 、 草 原 、 草 馈 、 训 漠 \ 沼 
泽 沙 生 植 被 ,栽培 植被 及 高 山 稀 跑 植被 ,植被 类 型 从 东南 问 西 北 , 依 次 为 森林 植被 、 和 森林 草原 植被 温 性 草原 
植被 荒漠 半 荒 漠 植 被 ” 。 


2 研究 方法 


2.1 资料 来 源 

本 文 所 用 数据 包括 黄土 高 原 2001 一 2014 年 MODIS-NDVI MODIS-LST 以 及 生态 区 划 和 植被 类 型 数据 。 
MODIS-NDVI 来 自 LP DAAC 中 的 MODIS 产品 的 MOD13A2 数据 集 ,空间 分 辩 率 为 1000m ,时 间 分 辩 率 为 16d， 
时 间 市 点 分 别 为 1、17v33、…、353; 利 用 单位 转换 系数 0.0001 将 其 转 为 0 一 1 之 间 ; MODIS-LST 来 源 于 
MOD11A2 ,空间 分 辨 率 为 1000m ,时 间 分 辨识 为 8d ,利用 单位 转换 系数 0.02 将 其 转 为 开尔文 单位 (Kelvin ) ,由 
于 LST 存在 大 面积 的 无 效 数 据 , 所 以 采用 Cressman 客观 分 析 法 与 非 对 称 高 斯 函数 拟 合 -2 ,将 8dLST 数据 进 
行 重 构 , 为 了 和 与 NDVI 数据 时 间 分 状 率 统一 ,采用 平均 值 合 成 法 将 其 合成 为 16d 数据 。 生 态 区 划 数 据 来 源 于 
中 国生 态 系统 评 佑 与 生态 安全 格局 数据 库 ( www.ecosystem.csdb.cn/index.jsp ) ,该 数据 集 是 在 生态 环境 问题 、 
生态 系统 敏感 性 .生态 系统 服务 功能 重要 性 的 基础 上 ,将 一 系列 相同 比例 尺 的 评价 图 ,采用 空间 芋 置 法 相关 
分 析 法 ~ 专家 集成 等 方法 , 按 生 态 功能 区 划 的 等 级 体系 ,通过 自 上 而 下 划分 方法 进行 生态 功能 区 划 方 案 划 分 ， 
生态 区 划 主 要 用 于 分 析 植 被 覆盖 与 土壤 湿度 年 内 关系 的 空间 分 布 以 及 植被 群落 特征 ,属于 植被 的 宏观 特征 。 
植被 类 型 数据 来 源 于 国家 科技 基础 条 件 平 台 一 国家 地 球 系 统 科 学 数据 共享 服务 平台 (http://www.geodata. 
cn) ,由 “七 五 ”重点 科技 攻关 项 目 (74-04-03-02) 成 果 系 列 图 扫描 后 数字 化 而 形成 500m 空间 分 辨 紊 的 栅 格 数 
据 集 ,主要 用 于 分 析 植 被 覆盖 与 土壤 湿度 年 内 关系 的 植被 类 型 特征 ,属于 植被 的 中 观 特 征 。 
2.2 土壤 湿度 估算 

Sandholt 等 指出 地 表 温 度 ( Land Surface Temperature，LST) 和 植被 覆盖 ( Normalized Difference Vegetation 
Index，NDVI) 的 梯形 空间 可 以 反映 土壤 水 分 高 低 ( 图 2) ,所 以 本 文采 用 TVDI 来 表征 土壤 湿度 ” ,TVDI 基于 
植被 指数 /地 表 温 度 的 梯形 特征 空间 计算 而 来 ,原理 如 下 : 
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Ts — Ts,,, 
TVDI = (1) 
Ts — Ts,,, 
其 中 ,7; 为 每 个 像 元 的 LST; Ts 为 对 应 像 元 湿 边 
的 LST, Ts 为 对 应 像 元 湿 边 的 LST( 图 2)。 Tamax 


干 边 和 湿 边 的 拟 合 是 计算 TVDI 的 关键 环节 。 首 
先 , 将 NDVI 的 值 域 0 一 1 划分 为 100 个 节点 ,节点 之 间 
相差 0.01; 然 后 ,以 每 个 节点 为 中 心 左右 分 别 扩 展 0.005 。 扫 
个 单位 形成 一 个 区 间 , 并 在 区 间 内 查询 LST 的 最 大 值 
和 最 小 值 ;最 后 ,分 别 以 节点 值 序列 为 自 变 量 ,LST 的 最 


植被 冠 层 
(水 分 胁迫 ) 


大 值 和 最 小 值 序列 为 因 变 量 进行 一 元 线性 回归 拟 合 ,从 “ 湿 一 (分 有 
而 得 到 干 边 和 湿 边 ,再 根据 式 (1) 计 算出 每 个 像 元 的 0 0.1 
TVDI。 本 文中 干 湿 边 拟 合 均 采用 Theil-Sen 趋势 法 ,以 NDVI 

此 剔除 数据 中 存在 的 异 浓 值 ,使 得 拟 合 的 干 湿 边 方程 更 


Ee A 图 2 ”NDVI 与 LST 构建 TVDI 示意 图 
为 稳健 。 0 值 域 在 < 之 加 2 IVDI 越 大 : 土壤 湿度 Fig.2 The TVDI configured by NDVI and LST 


越 低 ,TVDI 越 小 ,土壤 湿度 越 高 。 本 文采 用 1-TVDI 来 
表示 土壤 湿度 的 大 小 , 值 越 小 ,土壤 湿度 越 低 , 值 越 大 ,土壤 湿度 越 高 。 
根据 上 述 原 理 利用 16 d LST 和 NDVI 数据 计算 得 到 16 d TVBIL 数 据 , 值得 注意 的 是 ,TVDI 的 计算 相当 于 
是 对 LST 进行 最 大 限度 标准 化 ,标准 化 的 上 下 限 不 是 一 个 定 值 ,而 是 随 着 NBVL 的 不 同 发 生变 化 ,所 以 TVDI 
与 NDVI 并 非 同 源 数据 ,可 以 反映 出 土壤 湿度 与 植被 覆盖 之 间 的 关系 。 
2.3 ”数据 平滑 与 增 采 样 
计算 得 到 每 16d NDVI 和 1-TVDI 数据 的 多 年 均值 栅 格 序列 ,作为 年 内 NDVI 和 1-TVDI 时 间 序 列 数据 ,并 
对 其 采用 Whittaker 平滑 器 进行 平滑 一 ,原理 如 下 : 
f= ae) + 42 (Mf) (2) 
式 中 :y 为 原始 序列 ,/ 为 平 请 序列 ; 之 (xi = 了) 为 了 的 保 真 度 , A 为 正则 参数 , A” 为 阶 微分 算 子 。 为 了 方 
便 运 算 ,采用 线性 代数 计算 , 则 上 式 变 为 : 
f= gmidf[ (y Ue (3) 
FS +ADD) YY (4) 
式 中 :7 为 单位 矩阵 ;有 为 阶 差分 就 很 ,用 以 替代 4 。 记 工 =7+ AD:D ,L 越 大 (L=0,1,2,3),f 则 越 平滑 ， 
本 文 为 了 尺 可 能 保证 原 有 数据 变化 特征 ,数据 平滑 程度 采用 L=0。 
采用 德里 时 插值 法 平 请 数据 进行 增 采 样 。 该 方法 基于 快速 傅 里 叶 变 换 ( Fast Fourier Transform,，FFT) 将 
原始 16d 时 间 序 列 数据 变换 到 频 域 上 ,在 频 域 序列 中 扩充 采样 点 ( 补 零 ) ,然后 再 反 傅 立 叶 变换 , 增 采样 成 天 
时 间 序 列 (1 一 353d)。 
2.4， 时 请 互相 关 分 析 
本 文 基于 增 采 样 得 到 的 植被 覆盖 和 土壤 湿度 天 数据 ,根据 时 沛 互相 关系 数 原 理 计算 两 者 年 内 最 值 相关 系 
数 和 沛 后 时 间 ,计算 如 下 : 
(1) 首 先 ,计算 不 同时 滞 下 植被 覆盖 与 土壤 湿度 之 间 的 相关 系数 。 


> (4 — xi) (yi — Yin) 
ri(%,yY) = (9 ) 


ST -x,) ， SE yi ) 
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式 中 ,rm(x,y) 为 不 同时 请 天 下 植被 履 盖 与 土壤 湿度 的 相关 系数 序列 ;x; 与 刀 分 别 为 土壤 湿度 .植被 覆盖 序列 ;nm 
为 序列 长 度 六 为 滞后 时 间 (E=0,+1,+2,…) ,根据 经 验 其 值 小 于 等 于 n/4, 由 于 文中 是 对 天 数据 的 相关 分 析 ， 
所 以 n=353,k 的 最 大 值 取 90。 

(2) 计算 得 到 各 个 时 滞 下 的 最 值 互 相关 系数 以 及 其 对 应 时 滞 。 


n(ki)= max(ri(x,y)) (6) 
rk) = min(r(%,y)) (7) 
RS | | 
R=na ,KR =nal | |= "| | (8) 


R =r,,K=k,|r, | < [六 | 
式 中 ,为 时 滞 三 下 的 植被 覆盖 与 土壤 湿度 的 最 大 相关 系数 ;为 时 法 已 下 的 植被 覆盖 与 二 了 湿度 的 最 小 相 
关系 数 ;与 扩 分 别 为 植被 覆盖 与 土壤 湿度 的 最 值 互 相关 系数 及 其 时 灌 。 
K>0 表示 年 内 土壤 湿度 变化 影响 植被 覆盖 变化 ;K<0 表示 年 内 植被 覆盖 变化 影响 土壤 湿度 灰 化 ;K=0 表 
示 年 内 土壤 湿度 与 植被 覆盖 变化 同步 。 


3 结果 与 分 析 
3.1 黄土 高 原 植 被 覆盖 与 土壤 湿度 年 内 时 间 变 化 特征 0.6 0.50 
e 植被 覆盖 o 十 壤 温 度 | 


黄土 高 原 植被 绑 盖 和 土壤 湿度 年 内 季节 变化 ,如 图 
3 所 示 。 植 被 覆盖 年 内 呈 增 加 一 减 小 变化 。 最 低 值 出 
现在 第 17 天 (1 月 中 旬 ) ,之 后 逐渐 上 升 ,经 过 NDVI 的 
一 阶 差 分 可 知 植被 在 第 113 天 (4 月 下 旬 ) 左 右 生 长 发 
育 最 为 旺盛 ,植被 覆盖 变化 率 为 0.045/ 月 ;到 第 225- 天 

1 33 65 97 129 161 193 225 257 289321 353 

(8 月 中 旬 ) 植 被 覆盖 出 现 最 高 值 0.497 ,之 后 植被 开始 时 间 Time/d 
衰退 ,第 273 天 (10 月 上 名) 植被 衰退 最 为 明 显 ,植被 复 图 3 黄土 高 原 植被 覆盖 与 土壤 湿度 的 年 内 季节 变化 特征 
壮 变 化 率 为 -0.06L 月 。 Fig.3 The seasonal variations of vegetation cover and soil 

土壤 湿度 年 内 呈 减 小 一 增 如 他 小 入 六 增 放 变化 。 moisturein the Loess Plateau 
土壤 湿度 最 高 值 出 现在 一 年 内 的 1 月 上 旬 , 黄 土 高 原 处 
于 封冻 期 ,1 月 上 旬 至 4 月 下 旬 降 水 量 较 小 ,植被 逐渐 生长 发 育 并 达到 旺盛 期 ,植被 耗 水 使 得 土壤 湿度 逐渐 减 
小 至 最 低 值 ;4 月 下 旬 至 8 月 中 旬 降 水 量 增 加 ,植被 生长 至 最 大 覆盖 量 , 土壤 降水 补给 大 于 植被 生长 耗 水 , 植 
被 覆盖 与 土壤 湿度 均 不 断 增 大 ;8 月 中 旬 至 10 月 中 旬 ,植被 开始 退化 ,该 时 段 黄土 高 原 降 水 仍然 较 多 ,土壤 降 
水 补给 与 植被 生理 需 水 趋 于 平衡 ,土壤 湿度 平稳 变化 ,但 是 表现 为 微弱 的 减 小 变化 ;10 月 中 旬 至 年 末 , 植 被 进 
和 人 休眠 期 ,植被 生理 需 水 较 少 ,土壤 湿度 在 此 期 间 逐 渐 增 加 。 

可 以 看 出 ,植被 生长 对 土壤 湿度 变化 存在 明显 的 时 浸 效 应 ,植被 生长 以 及 植被 物候 与 土壤 湿度 变化 密切 
相关 ,植被 生长 季 始 期 与 土壤 湿度 谷 值 时 刻 吻 合 , 植 被 覆盖 峰值 时 刻 、 生 长 季 末 期 分 别 与 土壤 湿度 平稳 期 的 始 
期 与 未 期 一 致 ,反映 出 植被 不 同 生长 期 需 水 情况 ,以 及 植被 覆盖 变化 与 土壤 湿度 变化 之 间 的 相互 作用 。 

本 文 逐 个 顶 格 计算 得 到 2001 一 2014 年 植被 覆盖 和 土壤 湿度 的 多 年 16d 均值 数据 ,根据 时 清 互 相关 原理 ， 
提取 两 者 逐个 栅 格 的 年 内 时 间 序 列 ,经 过 Whittaker Smooth 平滑 和 健 里 叶 增 采样 将 数据 人 处理 为 天 数据 ,并 计算 
其 不 同济 后 时 间 对 应 的 相关 系数 ,根据 相关 系数 筛选 出 最 高 相关 程度 及 其 对 应 时 消 ( 图 4) ,来 解释 黄土 高 原 
年 内 植被 覆盖 与 土壤 湿度 的 相互 关系 。 

3.2 ”黄土 高 原 植被 覆盖 与 土壤 湿度 时 消 互 相关 空间 分 布 特征 

黄土 高 原 植被 覆盖 与 土壤 湿度 的 最 值 互相 关系 数 和 涡 后 时 间 空 间 上 分 异 明显 (图 4)。 最 值 互相 关系 数 

与 时 沛 空间 上 呈 明 显 的 负 相 关 性 ,空间 相关 系数 为 -0.84, 即 最 值 互 相关 系数 为 正 ,时 涉 为 负 , 反 之 则 正 。 
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黄土 高 原 农 业 与 草原 生态 区 北部 的 晋 中 一 陕 北 一 陇 东 高 原 丘 陵 半 干旱 草原 地 佛 .燕山 一 太行 山 山 地 落叶 
阔 叶 林 地 带 ,植被 覆盖 与 土壤 湿度 呈 负 相关 ,相关 系数 小 于 -0.85 ,土壤 湿度 变化 影响 植被 覆盖 变化 ,植被 覆盖 
产生 62 一 90d 的 滞后 效应 ,植被 生长 对 土壤 湿度 响应 迟缓 ,该 地 域 面 积 占 总 面积 的 20.73% ; 其 次 ,宁夏 平原 、 
河套 平原 分 属 荒漠 草原 地 带 和 典型 草原 地 带 , 均 有 黄河 流 经 ,植被 覆盖 与 土壤 湿度 的 最 值 相关 系数 介 于 
-0.55 一 -0.75 之 间 ,宁夏 平原 植被 覆盖 滞后 时 间 在 62 一 90d 之 间 ,河套 平原 介 于 0 一 62d ,面积 总 计 25.03%。 

正 相 关 区 域 仍 然 集 中 于 工 -12 ,但 河流 中 上 游 流 域 范 围 内 正 相 关 较 为 明显 , 目 东 辐 西 主要 包括 黄河 治 河 、 
泾 河和 渭河 ,相关 系数 大 于 0.75 ,土壤 湿度 滞后 于 植被 覆盖 在 30 一 69d, 自 西向 东 、 自 南 向 北 滞后 时 间 逐 渐 增 
加 ;温水 谷地 及 其 南北 两 侧 的 高 寒 草 甸 草原 地 带 , 相 关系 数 达到 0.85 以 上 ,土壤 湿度 清 后 于 植被 覆盖 小 
于 29d。 

综 上 ,黄土 高 原 植被 覆盖 与 土壤 湿度 的 年 内 相关 关系 及 其 滞后 效应 ,在 不 同 生态 区 表现 出 明显 的 差异 性 ， 
反映 出 不 同 的 植被 生态 群落 的 年 内 生长 与 土壤 水 分 相互 影响 特征 差异 性 明显 ;同一 生态 区 内 哩 具有 统一 性 ， 
但 不 同 地 域 上 也 表现 出 差异 性 ,黄土 高 原生 态 环 境 的 多 样 性 造成 生态 群落 内 部 还 具有 差异 性 ; 因此 ,本 文选 
用 植被 类 型 数据 分 析 了 不 同 植被 年 内 生长 与 土壤 水 分 的 相互 作用 关系 及 时 间 滞 后 效应 。 


最 值 互相 关系 数 滞后 时 间 
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图 4 黄土 高 原 植被 覆盖 与 土壤 湿度 最 值 相关 系数 和 滞后 时 间 
Fig.4 The ‘extrenium correlation coefficient and lag time of vegetation cover and soil moisture in the Loess Plateau 
生态 区 代码 I-10:; 蓝 山大 行 山 山 地 落叶 阔 叶 林 生 态 区 ,I-11: 汾 渭 分 地 农业 生态 区 ,1- 12: 黄 土 高 原 农 业 与 草原 生态 区 ,1-15: 秦 巴山 地 落叶 与 
常 绿 阔 叶 林 生态 区 , 工 - 五 内 蒙古 高 原 中 东部 典型 草原 生态 区 ,I-2: 内 蒙古 高 原 中 部 - 陇 中 苇 漠 草原 生态 区 ,I-3: 内 蒙古 高 原 中 部 草原 化 车 演 
生态 区 ,ITETL: 祁 连 山 森林 与 高 寒 草 原生 态 区 ,IIL-4: 江 河源 区 -甘南 高 寒 草 多 草原 生态 区 


3.3 “贡生 高 原 不 同 植被 类 型 植被 覆盖 与 土壤 湿度 时 浊 互 相关 分 析 

不 同 植被 类 型 的 植被 履 诲 与 土壤 湿度 相互 影响 有 所 差异 ,本 文 利用 黄土 蜗 原 不 同 植被 类 型 分 别 统计 最 值 
相关 系数 和 沛 后 时 间 的 平均 值 ,并 将 其 可 视 化 在 玫 现 图 中 ,如 图 5 所 示 。 

强 早生 好 类 草 ( 归 为 工 类 植被 ) 植 被 履 善 与 土壤 湿度 变化 呈 负 相关 ,相关 系数 为 -0.583 , 浏 后 时 间 平 均 为 
72d, 其 次 分 别 为 沼 渗 、 温 性 针 叶 林 落叶 阔 叶 林 与 温 性 落叶 灌 从 ( 工 类 植被 ) , 负 相 关 性 依次 减弱 ,植被 履 凑 产 
生 的 沸 后 效应 依次 减 小 至 11d( 图 5) ,说 明 这 5 种 植被 生长 对 土壤 湿度 变化 的 啊 应 程度 依次 减弱 ,但 啊 应 速度 
在 依次 加 快 , 主 要 是 由 于 土壤 湿度 年 内 变化 不 明显 引起 的 ;从 相关 系数 来 看 , 负 相 关 性 主要 表现 在 时 间 变 化 
上 :1 月 初 开 始 土壤 湿度 降低 ,到 3 月 中 下 旬 降 到 最 低 值 , 植 被 覆盖 变化 不 明显 ,之 后 植被 覆盖 随 土壤 湿度 增 
大 而 增 大 ,经 过 天 ( 沛 后 时 间 ) 后 植被 覆盖 达到 最 大 值 ,由 于 植被 耗 水 土壤 湿度 变化 平稳 ,在 此 期 间 及 其 后 
植被 覆盖 影响 土壤 湿度 变化 ,植被 覆盖 的 减 小 引起 土壤 湿度 增加 ,但 其 产生 的 负 相 关 歼 应 小 于 上 半年 土壤 湿 
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图 5 黄土 高 原 不 同 植被 类 型 植被 覆盖 与 土壤 湿度 最 值 相关 系数 和 滞后 时 间 

Fig.S The extremum correlation coefficient and lag time of vegetation cover and soil moisture of different vegetation types in the 
Loess Plateau 

名 称 颜 色相 同 的 植被 ,其 变化 特征 相似 ,被 归 为 同 二 类 ;其 中 工 类 植被 :红色 ; 开 类 植被 :黑色 ; 亚 类 植被 :绿色 ;区 类 植被 : 蓝 色 ; 无 植被 : 灰 
色 ; 特 殊 植被 竹林 :粉色 。 植 被 类 型 代码 :0 居民 起 冲 宕 温 性 针 叶 林 ;2 温 性 针 叶 林 ;3 落叶 阔 叶 林 ;4 杂 木 林 ;5 竹林 ;6 高 寒 落 叶 灌 从 ;7 高 寒 
常 绿 灌 从 ;8 温 性 落叶 灌 从 ;9 盐 地 落叶 灌 从 $10 草丛 ;11 草 甸 草原 ;12 典型 草原 ;13 荒漠 化 草原 ;14 高 寒 草原 ;15( 亚 ) 高 山 草 和 匈 ;16 河 、 湖 滩 
地 草 旬 ;17 低地 盐 生 草 多 ;18 强 旱 生 小 灌木 ` 半 灌木 ,灌木 ;19 强 旱 生 杂 类 草 ;20 沼泽 植被 ;21 沙 生 植被 ;22 草本 栽培 植被 ;23 木 本 栽培 植 
被 ;24 沙 地 ;25 高 山 稀 玻 植被 裸 岩 ;26 河滩 ;27 水 系 


度 影响 植被 覆盖 的 负 相 关 效 应 ,因此 总 体 上 表现 为 植被 覆盖 受 土壤 湿度 影响 (图 6) ;从 植被 覆盖 大 小 来 看 , 低 
履 盖 植被 生长 对 土壤 湿度 变化 的 啊 应 程度 较 高 , 啊 应 速度 较 慢 ,反之 高 瞻 盖 植 被 啊 应 程度 较 低 ,但 啊 应 速度 较 
快 ;从 植被 物候 来 看 , 随 着 植被 生长 季 长 度 增 加 ,植被 对 土壤 湿度 的 啊 应 程度 减 小 ,效应 速度 加 快 。 

低地 盐 生 草 旬 草丛 河滩、 寒 温 性 针 叶 林 、 河 一 湖 滩地 草 匈 、 木 本 栽培 植被 . 草 乞 草原、 盐 地 落叶 灌 从 与 杂 
木林 ( 下 类 植被 ) 植 被 禾 盖 与 土壤 湿度 的 负 相 关 性 依次 减弱 ,草本 栽培 植被 与 典型 草原 ( [[ 类 植被 ) 呈 正 相 关 ， 
但 其 平均 相关 系数 均 小 于 0.25 , 济 后 时 间 均 在 0 时 沛 附近 上 下 波动 (图 5) ,反映 出 植被 履 盖 与 土壤 湿度 同步 
变化 (年 内 变化 子 图 较 多 , 略 )。 

强 旱 生 小 灌木 一 半 灌 木 一 灌木 、 亮 江 化 草原 与 沙 生 植被 (于 类 植被 ) 均 值 相 关系 数 介 于 0.15 一 0.35 之 间 ， 
沛 后 时 间 小 于 20d(K<0) ,反映 出 其 土壤 湿度 变化 沾 后 于 植被 覆盖 变化 , 即 植被 覆盖 变化 影响 土壤 湿度 发 生 
变化 ,影响 的 强度 依次 增加 ,但 土壤 湿度 对 植被 履 盖 的 啊 应 速度 依次 变 慢 ( 图 5) ;从 相关 系数 来 看 , 正 相 关 性 
特征 表现 在 3 月 下 名 开始 土壤 湿度 随 植 被 覆盖 的 增加 而 增加 (图 7) ;从 植被 履 盖 大 小 来 看 ,三 者 均 属于 低 窗 
盖 植 被 ,植被 窗 六 年 内 变化 较为 平 绥 但 其 最 值 依次 增 大 ,植被 儿 盖 与 土壤 湿度 变化 特征 依次 变 得 明显 ,从 而 表 
现 为 植被 覆盖 对 土壤 湿度 的 影响 依次 增强 。 
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图 6 I 类 植被 NDVI 与 1-TVDI 的 年 内 季节 变化 
Fig.6 The seasonal variations of NDVI and 1-TVDI of vegetation I in the Loess Plateau 
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图 7 ”了 H 类 植被 NDVI 与 1-TVDI 的 年 内 季节 变化 
Fig.7 The seasonal variations of NDVI and 1-TVDI of vegetation LI in the Loess Plateau 


高 山 稀 巩 植被 一 裸 罕 、( 亚 ) 高 山 草 多 .高 寄 稼 绿 洪 从 .高寒 落叶 灌 从 与 高 寒 草 原 (IV 类 植被 ) 均值 相关 系 
数 依次 增 大 ,并 介 于 0.6 一 0.75 之 间 ,滞后 时 间 介 于 0 一 30d(CK<0)( 图 5) ,说 明 高 山 / 高 寒 植被 生长 对 土壤 湿度 


变化 的 


啊 程 度 最 强 ,而 且 啊 应 速度 较 快 ;从 相关 系数 来 看 , 正 相 关 性 依次 增加 , 主要 是 由 于 土壤 湿度 年 内 变 


化 呈现 出 增加 一 减 小 分 布 ,与 植被 覆盖 变化 基本 吻合 , 且 土 壤 湿度 变化 特征 越 来 越 明 显 , 表 现在 时 间 变 化 上 ,1 
月 初 土壤 湿度 随 植 被 履 凑 逐渐 增加 ,到 6 月 中 下 旬 植 被 达到 最 大 罗 善 ,植被 生理 需 水 较 高 致使 土壤 湿度 减 小 ， 
植被 事 退 后 土壤 湿度 有 所 回升 ,到 8 月 上 中 旬 土 壤 湿 度 随 植被 窗 盖 减 小 而 减 小 (图 8) ;从 植被 覆盖 大 小 来 看 ， 
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植被 生长 对 土壤 湿度 变化 的 影响 程度 ,及 土壤 湿度 的 啊 应 速度 ,与 植被 覆 将 大 小 关系 不 明显 ,土壤 湿度 的 啊 应 
速度 由 快 到 慢 依 次 为 :高 寒 章 原 >( 亚 ) 高 山 草 向 > 蜗 宕 落叶 灌 从 > 局 守 常 绿 灌 从 > 高 山 稀 玖 植被 一 昼 宕 。 

值得 注意 的 是 ,人 竹林 生长 需 水 条 件 高 且 敏 感 ,致使 植被 覆盖 与 土壤 湿度 均值 相关 系数 高 达 0.55 , 沛 后 时 间 
平均 为 20d( 图 5) ,植被 覆盖 变化 沛 后 于 土壤 湿度 变化 。 
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图 8 TV 类 植被 NDVI 与 1-TVDI 的 年 内 季节 变化 
Fig.8 The seasonal variations of NDVI and 1-TVDI of vegetation IV in the Loess Plateau 
植被 类 型 代码 ;6 高 寒 落叶 灌 从 ;7 高 寒 常 绿 灌 从 ;14 高 宕 草原 ;15( 亚 ) 高 山 草包 ;25 高 山 稀 焉 植被 . 裸 宕 


4 讨论 


土壤 湿度 年 内 季节 变化 特征 与 张蕾 等 的 西北 地 区 土壤 湿度 变化 特征 ”高 度 吻 合 ,差异 性 主要 表现 在 土 
壤 湿 度 变 化 特征 点 所 对 应 的 时 间 点 不 同 。 上 土壤 湿度 最 高 值 出 现在 一 年 内 的 1 月 上 旬 ,4 月 下 名 土壤 湿度 降 至 
最 低 值 ,而 张 蔷 等 分 析 土 壤 湿 度 12 月 中 旬 至 1 月 较 高 ,5 月 下 旬 至 6 月 上 旬 降 至 最 小 ,本 文 土壤 湿度 最 低 
值 对 应 时 间 较 张 雷 提前 了 1 个 月 以 上 。 原 因 主 要 是 数据 的 时 间 分 辩 率 不 同 引 起 的 ,采用 16d 数据 进行 分 析 ， 
而 张 一 采 用 旬 值 数据 ,分 析 中 体现 了 植被 生长 耗 水 引起 土壤 湿度 变化 ,虽然 其 时 间 分 辩 率 较 高 ,但 缺乏 植被 柳 
蓄 数 据 作为 支撑 ,土壤 湿度 特征 时 间 市 点 与 实际 植被 生长 情况 不 符 。 分 析 了 植被 禾 盖 年 内 季 市 变化 特征 , 植 
被 生长 季 始 期 为 4 月 下 旬 , 与 土壤 湿度 最 低 值 时 间 点 契合 度 很 好 , 韦 振 锋 等 分 析 得 出 陕 甘 宁 黄土 高 原 大 多 数 
植被 生长 季 始 期 为 3 一 4 月” ,李强 等 认为 生长 季 始 期 为 3 月 下 旬 至 5 月 下 旬 ” , 谢 宝 妮 等 认为 生长 季 始 期 
为 了 月 卡 旬 至 ,5 月 下 名 ,证 明 TVDI 用 来 表示 土壤 湿度 变化 特征 具有 很 高 的 可 信 性 。 其 次 ,本 文 在 利用 时 
沛 互相 关 方 法 之 前 , 随 年 内 植被 覆盖 和 土壤 湿度 数据 进行 了 Whittaker 平滑 和 傅 里 叶 插 值 增 采样 ,一 定 程度 上 
别 除 也 数据 噪声 ,提高 了 信 品 比 。 

从 不 同 植被 类 型 来 看 ,植被 生长 与 土壤 水 分 的 年 内 变化 具有 和 密切 的 关系 ,土壤 水 分 的 增加 可 以 促进 植被 
生长 ,而 植被 生长 又 消耗 土壤 水 分 ,从 而 曾 弱 了 土壤 水 分 增加 速度 。 对 于 非 高 山 、 非 高 晨 植被 来 说 ,上 半年 土 
二 湿度 的 谷 值 是 返青 期 植被 快速 生长 引起 的 ,下 半年 土壤 湿度 的 谷 值 是 植被 从 戊 盛 期 到 枯黄 期 为 了 维持 基本 
的 生理 机 能 ,摄取 的 土壤 水 分 大 于 雨水 等 补给 造成 的 ;对 于 高 山 高 寒 植被 来 说 ,植被 生长 需 水 较 少 ,植被 覆 装 
与 土壤 湿度 增加 同步 ,返青 期 植被 快速 生长 时 ,土壤 水 分 大 量 耗 减 ,过 了 植被 戊 盛 期 ,土壤 水 分 因为 各 种 补给 
以 及 植被 生理 需 水 减 小 而 增加 ,年 末 补 给 减 小 土壤 水 分 也 随 之 减 小 ,从 而 形成 两 个 峰值 。 黄 土 高 原 土 壤 水 分 


http://www.ecologica.cn 


6 期 张 囊 ”等 :黄土 高 原 植被 覆盖 与 土壤 湿度 的 时 沛 关联 及 时 空 特 征 分 析 2137 


年 内 变化 的 峰值 和 谷 值 与 植被 物候 期 非常 吻合 ,因此 土壤 水 分 的 年 内 变化 可 作为 植被 物候 特征 提取 的 一 种 重 
要 依据 。 

由 基于 沛 后 时 间 的 最 值 互 相关 系数 计算 原理 可 知 , 该 方法 可 以 在 一 定 程 度 上 反映 植被 窗 盖 和 土壤 湿度 的 
年 内 相互 影响 机 制 。 但 是 当 最 大 和 最 小 相关 系数 差异 性 较 小 ,尤其 相关 系数 分 别 为 正 值 和 人 负 值 时 ,得 到 的 最 
值 相 关系 数 不 能 客观 反映 植被 窗 盖 和 土壤 湿度 关系 , 祁 后 时 间 难 以 区 分 植被 履 新 影响 土壤 湿度 变化 还 是 土壤 
湿度 影响 植被 禾 盖 变化 ,因此 存在 误 判 的 情况 。 为 外 ,植被 物候 特征 与 土壤 水 分 变化 密切 相关 ,植被 生理 需 水 
状况 可 能 是 两 者 关系 的 契合 点 ,以 后 应 在 考虑 植被 不 同 生理 阶段 消耗 土壤 水 分 的 基础 上 分 析 植 和 被覆 冀 与 土壤 
湿度 的 相互 关系 。 


5 结论 


(1) 黄土 高 原 植 被 后 长 对 土壤 湿度 变化 存在 明显 的 时 淖 效 应 ,植被 覆盖 与 土壤 湿度 之 间 的 相互 作用 渗透 
于 植被 生长 的 整个 过 程 ,植被 生长 以 及 植被 物候 与 土壤 湿度 变化 密切 相关 ,反映 出 不 同 植 被 生长 期 需 水 情况 。 

(2) 黄 土 高 原 植被 窗 六 与 土壤 湿度 之 间 的 最 值 互 相关 系数 与 时 沛 的 空间 分 布 旺 负 相关 。 黄 土 高原 土 壤 
湿度 对 其 植被 窗 六 变 化 影响 强烈 区 ,主要 分 布 在 农业 与 草原 生态 区 北部 的 丘陵 沟 蜜 区 ,但 植被 生长 对 土壤 湿 
度 啊 应 迟缓 ;植被 禾 冀 对 土壤 湿度 变化 影响 强烈 区 ,分 布 于 西部 的 高 宕 地 带 , 啊 应 时 间 相 对 最 短 , 其 次 为 农业 
与 草原 生态 区 内 河流 中 上 洲 流 域 , 且 土 壤 湿 度 对 植被 覆盖 的 啊 应 较 快 ; 目 东 癌 西 月 北 辐 南 啊 应 逐渐 加 快 。 

(3) 强 早生 杂 类 草 .沼泽 . 温 性 针 叶 林 、 沙 叶 阁 叶 林 与 温 性 落叶 淮 俱 的 植被 履 状 与 土壤 湿度 呈 较 强 负 相 关 
性 ,植被 生长 对 土壤 湿度 变化 的 啊 应 程度 依次 减弱 ;响应 速度 依次 加 快 。 低 地 盐 生 草 旬 .草丛 河滩、 寄 温 性 针 
叶 林 河 一 湖 滩 地 草 甸 、 木 本 栽培 植被 . 草 旬 草原 . 盐 地 洛 叶 浇 傣 条 木林 .草本 栽培 植被 与 典型 草原 的 植被 履 
闭 与 土壤 湿度 的 相关 性 很 弱 , 浏 后 时 间接 近 于 0d ,反映 出 植被 复 盖 与 土壤 湿度 同步 变化 。 土 壤 湿度 变化 济 后 
于 强 旱 生 小 灌木 一 半 灌 木 一 灌木 . 充 演 化 草原 与 沙 生 植 被 的 植被 禾 盖 变化 ,植被 覆盖 正 问 影 响 土壤 湿度 变化 
的 强度 依次 增加 ,但 时 间 效 应 依次 延长 。 高 山 / 高 宕 植 伞 生 长 对 土壤 湿度 变化 的 正 问 影响 程度 最 强 ,而 且 啊 应 
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